Pseudo-Boundaries in Discontinuous 2-Dimensional Maps 
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It is known that Kolmogorov-Arnold-Moser boundaries appear in sufficiently smooth 2- 
dimensional area- preserving maps. When such boundaries are destroyed, they become pseudo- 
boundaries. We show that pseudo-boundaries can also be found in discontinuous maps. The origin 
of these pseudo-boundaries are groups of chains of islands which separate parts of the phase space 
and need to be crossed in order to move between the different sub-spaces. Trajectories, however, do 
not easily cross these chains, but tend to propagate along them. This type of behavior is demon- 
strated using a "generalized" Fermi map. 



The most obvious way to set particles into motion is to 
subject them to an external field. Yet, particle currents 
can be also induced by zero averaged time-dependent 
forces, provided an asymmetric potential is applied in the 
system. Recent studies of systems consisting of Brownian 
particles in periodic, locally asymmetric, potentials have 
shown that motion may be induced by periodic forcing, 
non-random noises or due to an oscillatory change in the 
potential profile Q. Such rectification processes are in- 
teresting because they can lead to the generation of new 
microscopic pumping techniques. They may also provide 
some insight into the mechanism of protein motors |2j. 
Another example for a net directed motion of particles 
in periodic (both in time and space) potentials is found 
in a model describing the motion of a classical particle 
within an infinite set of equally spaced potential barriers, 
whose position and height oscillate periodically M. At 
each collision, the particle either crosses or reflects from 
the barrier, depending whether or not its kinetic energy 
is larger than the potential height. The special feature of 
this system (unlike these in Refs. Q) is its chaotic deter- 
ministic nature. Chaotic dynamics are frequently studied 
by investigating the corresponding trajectories in the sys- 
tem's phase space. Our system is one-dimensional with 
a time dependent Hamiltonian, hence it's phase space 
is four-dimensional. The Poincare method can be ap- 
plied to reduce this 4-dimensional phase space into a 2- 
dimensional one by taking a surface lying in the phase 
space and studying the sets of intersections of different 
trajectories with this surface. Successive crossing points 
are computed using a 2-dimensional area-preserving map 
of the surface, the reduced phase space, onto itself. In or- 
der to inquire the features of this map, we constructed a 
picture of it's 2-dimensional phase space by plotting the 
sets of points (which are the trajectories in the reduced 
phase space) starting at different initial conditions J|||. 

In this paper we report and qualitatively explain an in- 
teresting phenomenon observed while studying this map. 
Pictures of the phase space, constructed for different val- 
ues of the map's parameters, displayed that the phase 
space consists of a stochastic sea, which a typical trajec- 
tory covers ergodically and some embedded islands, in 
which the motion is mainly regular over closed, elliptic- 
like, curves. This structure is quite common to many 
area-preserving maps. We did not observe Kolmogorov- 



Arnold-Moser (KAM) curves in our phase space. These 
curves of regular motion, if they exist, bound different 
parts of the stochastic sea, that is, prevent trajectories 
from leaving or entering them. According to the KAM 
theorem Ij-g], KAM curves do not appear in discontin- 
uous maps as is ours H. The origin of discontinuity in 
our map can be easily understood: Points in the phase 
space are divided into two classes each corresponding to 
either crossing or reflection from barriers. For each class 
the map is represented by one of its two possible "sub- 
maps" . Discontinuity appears on the boundaries between 
regions containing points of different classes (see Fig. [j]). 

Since KAM curves do not exist, the stochastic sea 
is a single connected unit with only a few islands ex- 
cluded pi. The motion over a connected (invariant) 
component of the phase space is ergodic and therefore a 
trajectory starting at almost any point in the stochastic 
sea, will eventually reveal this entire region. Numerical 
observations have shown that the number of iterations 
needed to explore completely the stochastic sea varies, 
within more than 3 orders of magnitude, for different 
values of the map's parameters. In the cases when this 
number was considerably large, the stochastic sea seemed 
to be divided into smaller parts in which trajectories trav- 
eled for long time intervals before moving from the one 
to the other. The boundaries between the different parts 
are called pseudo (partial) boundaries. Their appearance 
does not violate the ergodicity hypothesis, it simply im- 
plies that it is valid only for time scales much larger than 
the characteristic escape time from each pseudo-bounded 
part of the stochastic sea. 

Each KAM boundary is destroyed when the map's con- 
trol parameter exceeds some critical value. At the crit- 
ical point an infinite number of gaps are formed in the 
torus and it turns into a Cantor set (hence called a can- 
torus) JnLO]. As we continue to change the control pa- 
rameter, the flux across the cantorus increases from its 
zero value at the critical point [[13]]. Thus, close to the 
critical point, the flux across a cantorus is low and can- 
torus serves as a pseudo-boundary. Cantori, like KAM 
tori, do not appear in discontinuous maps. We show in 
this paper that pseudo-boundaries are found in discon- 
tinuous maps due to other reasons. They are formed by 
chains of islands surrounding different areas of the phase 
space. In order to leave these areas, trajectories need to 



diffuse across the chains, which usually appear in groups. 
However, the motion in the vicinity of the chains, tends 
to be along the chains, in a way that resembles the reg- 
ular motion over the chains. 

In smooth maps, cantori and chains of islands are 
deeply related. For each torus, cantorus or chain of is- 
lands, there is a corresponding winding number (for def- 
inition of winding numbers see, for example, |14[| ). Tori 
and cantori have irrational winding numbers while chains 
of islands have rational ones. Green jl5| has shown that 
close to any torus or cantorus, there are infinitely many 
chains of islands whose winding numbers are successive 
finite truncations of the continued fraction representing 
the irrational winding number of the torus or cantorus. 
It has also been shown [tl€3|Jl Ti l that the flux across a can- 
torus is the limit of the flux across its approximating 
chains of islands. This last result means that chains of 
islands which approximate cantori that serve as pseudo- 
boundaries (i.e., which the flux across them is low) are 
pseudo-boundaries themselves. For discontinuous maps, 
in which winding numbers are meaningless, the proper- 
ties of the chains should be studied differently, since they 
appear with no relation to cantori. A special emphasis 
should be given to the way in which trajectories behave 
when they are trapped between several adjacent chains. 

We consider a classical particle of unit mass, moving 
within a one-dimensional sequence of equidistant poten- 
tial barriers, each of an infinitesimal width and finite 
height, H(t). The barriers oscillate with the same fre- 
quency and in phase. Their velocity is given by v(t) = 
v b g(t), while H(t) = H [l + g(t)], where g(t) = sin(27rf). 
The particle moves freely between the barriers, occasion- 
ally colliding with them. At each impact it can either 
cross or be reflected from the barrier. It crosses the bar- 
rier if its kinetic energy in the reference frame of the 
barrier exceeds the height of the potentia l barrier at the 
moment of impact, i.e., \V n — v(t n )\ > y / '2H(t n ), where 
t n and V n are, respectively, the time of the n-th impact 
and the velocity of the particle before that impact. If the 
particle crosses the barrier, its velocity is not changed, 
i.e., V n +i = V n . If it is reflected, on the other hand, it 
acquires a new velocity: V n +\ = —V n + 2v(t n ). In order 
to calculate the moment of the next impact, one needs to 
consider the motion of both the particle and the barriers. 
We made the approximation that the distance traveled 
between two consecutive collisions is constant and equals 
to the distance between to barriers, L, as if the barri- 
ers do not change their position. We also assumed that 
in a case of reflection, the particle is reflected backward 
and therefore if V n + 1 has the same sign as V n , we re- 
place it with — V n +\. A similar approach was used in 
Rcf. |L8|] for the Fermi accelerator model (a ball bounc- 
ing between two walls [|9[). Due to the periodicity of 
g(t) only time modulo period (=1), which we denote as 
phase, is relevant, and two consecutive phases are related 
by t n+1 = t n +L/\V n+1 \ (mod 1). 

We chose the mean height of a barrier, An, to be the 
control parameter, while we set the other parameters 



L = 100 and v b — \. For the special case H = oo, 
barrier crossing is impossible and the map reduces to the 
Fermi map |g(J. For finite Hq, the map is discontinuous. 
Fig. [I] shows the discontinuity lines in the (V — t) phase 
space which separate the regions mapped by the differ- 
ent sub-maps of crossing and reflection (for Hq — 500). 
The map is area-preserving since both sub-maps are area- 
preserving and since their ranges do not overlap. For 



V > v m 



2V#o + v b and V < V m 



-2^/Ho + v b , 



points in the (V — t) phase space are mapped only by 
the crossing sub-map (the kinetic energy is sufficiently 
large to cross the barriers at any phase), hence the mo- 
tion in this part of the phase space is over the lines 
V = Vq = const. The stochastic sea lies between 
Vmin an d Vinax, including several islands. As already ex- 
plained, a picture of it is made by starting at almost any 
one of its points, iterating the map a sufficiently large 
number of times and plotting the resulting trajectories, 
{(t n ,V n );n = 0,1,2...}. 
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FIG. 1. Regions in the (t — V) phase space of points mapped by 
the sub-maps of crossing and reflection, for a barrier mean height 
Ho = 500. The map is discontinuous over the curves separating 
between the different regions. 

We performed this procedure for various values of Hq. 
For Hq values below 1000, the stochastic sea was com- 
pletely explored up to the resolution of the plots) after 
10 6 — 10 7 iterations. As expected, apart from some em- 
bedded chains of islands, the trajectory wandered in the 
entire region between Vmin and Vm ax - For H > 1000, the 
situation was quite different. Fig. [| depicts the velocities 
V n and phases t n of first 2 • 10 10 collisions between the 
particle and the barriers. The velocities appear to remain 
smaller than 63, while the stochastic sea ranges between 
Vni„ = -99.5 and V max = 100.5. 

A pseudo-boundary exists around \V\ ~ 63. A closer 
look of that area of the phase space reveals a few ad- 
jacent chains of islands, located one above the other, 
whose shapes are similar to the boundary of the area 
plotted in Fig. g. One of these chains is shown in Fig. ||. 
It is composed of two velocity branches, in the upper 
(V > 0) and the lower (V < 0) halves of the phase space, 



between which the motion alternates. The effect of this 
chain on the motion close to it is depicted in the inset 
of Fig. ||, showing an enlargement of a small segment of 
the chain. A narrow stochastic layer appears below the 
chain (above it, if the negative velocity branch is con- 
sidered). We found that trajectories can move over this 
narrow layer for thousands of iterations before leaving 
it, usually falling back into the stochastic sea. Two fea- 
tures resemble the behavior close to a KAM torus or a 
cantorus: (1) the motion is restricted to one side of the 
boundary (the chain), and (2) the motion is correlated 
to the boundary and only gradually deviates from it. As 
for a cantorus, if the chain is long and narrow with only 
small gaps between the islands, there are only small areas 
from which a trajectory can jump from one of its sides 
to the other. 
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FIG. 2. The subspace explored after N = 2 • 10 iterations of 
the map for a barrier mean height Ho = 2500. The points denote 
the velocity of the particle V n (y-axis), and the phase of the barrier, 
t n (x-axis), at the ra-th (n = 1, 2, . . . , N) collision. 

Fig. [^ demonstrates the effect of a single long chain 
with islands located close to each other. The stochastic 
motion near such a chain is highly correlated. It, roughly, 
follows the regular motion over the chain and only slowly 
deviates from it. But when a pseudo-boundary appears, 
we usually observe a group of several chains, lying closely 
one above the other, and not only a single chain. Trajec- 
tories need to diffuse across all these chains in order to 
move to another sub-space. However, due to the charac- 
ter of the motion near the chains, they tend to move along 
them. Suppose a trajectory enters between a few chains, 
close to one of them. It propagates along the chain, while 
gradually moving away. If there are a few nearby chains, 
than while it moves away from one chain it gets closer 
to another, along which it propagates for subsequent it- 
erations. This is an over-simplified picture since it is not 
always possible to relate the motion to a particular chain 
at each instance, yet, it qualitatively explains the charac- 
ter of the motion in this "band of chains" . The trajectory 
can be trapped between the chains for many iterations. 
If it succeeds in crossing all of them, it escapes. Usually, 
however, after spending some time wandering between 



the chains, it "falls back" into the stochastic sea. Fig. [| 
depicts six chains and two stochastic trajectories. Note 
that while one of the trajectories (open squares) travels 
between the chains labeled 1 to 5, the other trajectory 
(dots) is restricted to the area above the chain labeled 5. 
Indeed, since the flux across this chain is very small, com- 
pared to the flux across the other chains (probably be- 
cause of the fact that it consists of relatively many islands 
with small gaps between them), the chain labeled 5 (also 
shown in the inset of Fig. ||) serves as the natural limit 
for both sides of the pseudo-boundary. 
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FIG. 3. A chain which appears in phase space around \V\ ~ 63 
for a barrier mean height Ho = 2500. Its shape resembles the 
boundary of the area shown in Fig. El The inset is a magnification 
of a small part of the chain, showing few of its islands. A narrow 
stochastic layer is located just below the chain which serves as a 
pseudo-boundary. Trajectories that reaches the stochastic layer, 
may move over it for thousands of iterations before either falling 
back into the stochastic sea (the usual case) or crossing the chain 
(the rare scenario). 



In conclusion, we have shown that long and narrow 
chains of islands may serve as pseudo-boundaries in dis- 
continuous 2-dimensional maps. If a group of several ad- 
jacent chains surround a certain sub-space, they all need 
to be crossed in order to move to a different sub-space. 
However, in the vicinity of the chains, trajectories which 
stochastically wander between the islands tend to prop- 
agate along the chain, slowly moving from one chain to 
the other. Hence, the chains are not easily crossed and 
they form a pseudo-boundary. 
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FIG. 4. A few adjacent chains of islands found in the phase 
space for a barrier mean height Ho = 2500, forming together a 
pseudo-boundary. The numbers inside the islands label the differ- 
ent chains. The number of islands in the chains labeled 1 to 6, 
in the upper (lower) half of the phase space, are respectively: 27 
(33), 37 (47), 49 (59), 60 (72), 420 (466), 157 (164). The dots and 
the open squares belong to two stochastic trajectories wandering, 
respectively, above and below the chain labeled 5. (Only a small 
part of the positive velocity branches of the chains is shown here.) 
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